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Abstract 

A purely group-theoretical approach (for which the symmetric group plays a central 
role), based upon the use of properties of fractional-parentage coefficients and isoscalar 
factors, is developed for the derivation of the Coulomb energy averaged over the states, 
with a definite spin, arising from an atomic configuration n£^ . 
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1. Introduction 



In the framework of an alternative parametrisation for the theory of complex spectra 
[1-3], Kibler and Katriel [2] have conjectured a formula, viz, 

■(^ + 1)- 5(5 + 1)1 A (1) 

giving the average Coulomb energy Eg,y[n£^] for the states, with a fixed value of the spin 
S, which arise from an atomic configuration ni^ of equivalent electrons in a shell n£. 
In Eq. (1), the parameters S and A are defined through 

^EavM^] = E + A ^i^avN'] = S (2) 

corresponding to the singlet {S = So = 0) and triplet {S = So = I) states for the case 
N = 2. A formula equivalent to (1) and (2) has been derived independently by using 
the method of moments [4]. 

It is the aim of this short paper to give a straightforward proof of (1) by making 
use of fractional-parentage coefficients and of an unusual sum rule satisfied by isoscalar 
factors for a chain of compact groups closely connected to the permutation group. In- 
deed, we work out (in Section 2) an extension of formula (1), valid in the general case 
of a spin-independent two-body interaction, and then specialise it to the case of the 
Coulomb interaction. Some concluding remarks are given in Section 3. The relevant 
material (Racah's lemma and orthogonality-completeness property for Clebsch-Gordan 
coefficients or isoscalar factors) concerned with the group-theoretical approach of Section 
2 is relegated in an appendix. 



2. Average Energies 
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In the Russell- Saunders coupling (or LS'-coupling), the state vectors of the electronic 
configuration n£^ (for atoms and ions) may be written as 

* = \n£^[f]aLSMLMs). (3) 

Most of the symbols in (3) have their usual meaning. Let us simply precise that [/] 
stands for a Young pattern, with two columns (say of lengths (pi and (p2 with ipi > ip2), 
which characterises the orbital part of The spin part of may be described by the 
Young pattern [cp] = [(pi(p2], with two rows of lengths (pi and (p2, which turns out to be 
the transposed pattern [/] of [/]. The pattern [(p] is (unambiguously) connected with 
the number of electrons N and the total spin S via 

(pi + (p2 = N (pi-(p2 = 2-5. (4) 

Finally, the symbol a in (3) collectively denotes the remaining labels that are necessary 
for a single- valued enumeration of the (allowed) atomic state vectors of the configuration 

The vectors are expressed in a basis adapted to the following chain of groups 

U41+2 D U2£+l{-^ SO3 ^ SO2) ® U2{-^ SU2 ^Ui). (5) 

Each of the vectors (3) spans the antisymmetric irreducible representation class (IRC) 
{l"^} of the unitary group The Young patterns [/] and [p] characterise IRC's 

of the groups U2(.+\ and U2 for the orbital and spin parts of respectively. The 
quantum numbers L and S refer to IRC's of the subgroups SO3 and SU2 of U2£-\-i and 
U2, respectively. Finally, the projections and Ms indicate in turn IRC's of the 
subgroups SO2 and C/i of SO^ and SU2i respectively. 
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For a given configuration n£^, the total number of the vectors (3) allowed by the 
Pauli principle coincides with the dimension of the IRC {l^} of the group It is 

thus given by dim {l^} = C|g_|_2 in terms of binomial coefiicients. Among these C^^^2 
state vectors, we select the ones having a fixed value S of the spin. Obviously, the 
number of vectors with S fixed is nothing but the product Nf{2S + 1), where Nf is 
the dimension of the IRC [/] of U2e+i and 2S + 1 may be visualised as the dimension 
% of the IRC [ifi] of C/2. 

In the present paper, we are interested in the average energy 



^^+^EavK^] = ^ E E E {ni''[f]aLSMLMs\V\ne''[f]aLSMLMs) 

f aL Ml Ms 

(6) 

where F is a spin-independent two-body Hamiltonian 

TV 

v= with v,,^v,jin,). (7) 

j>i=i 

Equation (6) gives the average of V over the Nf{2S + 1) state vectors (3) having a 
fixed spin S. In Eq. (7), the sum on i and j is to be extended over the N{N — l)/2 
two-electron interactions Vij ; furthermore, Vij depends only on the distance between 
the electrons i and j. 

Invariance of V under the rotation group SO 3, and the spin group SU2 ensures that 
(6) can be rewritten as 



^^+^E^y[n£^] = "S2{2L + 1) {n£^aLS\V\n£^aLS). (8) 

J aL 

The next step is to introduce coefficients of fractional parentage in order to calculate 



the matrix elements of V in (8). This leads to [5] 




{raLS{\e 



-^a'L'S',fLoSo) 



{ne^LoSo\Vi2\ne^LoSo) 



Following Racah [6, 7] and Neudatchin and Smirnov [8], the two-particle coefficient of 



fractional parentage {\e^~^,e^) = {£^-^,£^1} e^)* in (9) can be developed as 



{£''aLS{\£''-^a'L'S'J^LoSo) = ^ 

(£^[/]aL{|£^-2[/V'L',£2[/,]Lo) (10) 
{s''[f]S{\s''-'[f']S',s'[fo]So) 



where s = 1/2. In Eq. (10), the IRC's [/], [/'] and [/o] refer to the group U2e+i- However, 
the symbols n/, n// and n/^ denote the dimensions of [/], [/'] and [/o] as considered 
as IRC's of the permutation groups Sn, Sn-2 and 5'2, respectively. From a group- 
theoretical viewpoint, the fractional-parentage coefficient {£^ aLS{\£^~^a' L' S' , £'^LoSo) 
is identical to the isoscalar factor {{l^-'^}a' L' S' + {l'^}LoSo\{l^}aLS) for the chain 
U4e+2 D SOs^SU2 [7]. Equation (10) then corresponds to the factorization (see Refs. [6] 
and [7]) 



{{l''-'}a'L'S' + {l'}LoSo\{l''}aLS) = ({l^-^}[/'][/'] + {l^}[/o][/o]|{l^}[/][/]) 



in terms of isoscalar factors for the chains D C^2^+i <8) U21+1 13 SO^ and 




([/>'L' + [/o]Lo|[/]aL) 



ms' + [h]So\[f]S) 



(11) 
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U2 D SU2. Indeed, the correspondence between (10) and (11) yields 



+ {l^}[/o][/o]|{l^}[/][/]) = 

V ^/ 

([/V'L' + [fo]Lo\[f]aL) = {e^[f]aL{\e^-'[f']a'L',e^[fo]Lo) (12) 

mS' + [/o]^o|[/» = {s'^[f]S{\s'^-'[f']S\s'[fo]So). 
In Eqs. (10) and (12), it is clear that n/^ = 1 for both IRC's [fo] = [2] and [/o] = [11] of 
82- In the case [fo] = [2] only even values Lq = 0,2, ■ ■ ■ ,2£ are admissible and Sq — 
while for [fo] = [11] we have Lq = 1, 3, • • • , 2£ — 1 and Sq = 1. Furthermore in view of 
(12), the spin part {s^[f]S{\s^~'^[f']S',s'^[fo]So) of the fractional-parentage coefficient 
(10) is trivial in the atomic case since the isoscalar factor {[f']S' + [fo]So\[f]S) for the 
chain U2 D SU2 is unity for the allowed values of S' and Sq. By combining Eqs. (8)-(12), 
we obtain 



2S+1 



^ aL a'L' S' Lq Sq 



^ (2L + 1) {n£''LoSo\Vi2\ne^LoSo) \{[f']a' L' + [fo]Lo\[f]aL)\\ 

Uf 

(13) 

An important step is now to effectuate the summations over aL and a'L' in (13). This 
may be easily done by using the orthogonality-completeness relation (22) (see appendix) 
applied to the (orbital) isoscalar factors for the chain U2e+i D SO^. We thus get 

S' Lq So (14) 

A([/]|[/'] ® [fo]) ^ (2Lo + 1) {ne^LoSo\V,2\nfLoSo). 
By introducing in (14) the average energy 

^^°+^E,,[ne^] = J2{2Lo + l)ine^LoSo\V^2\ne''LoSo) (15) 
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for the configuration nP, we finally arrive at 

= i Ar(Ar - 1) J] J] 

S' or [/'] So or [/o] (^26) 
^^o+^Ea^K^] A([/]|[/']®[/o])^. 

Uf 

Equation (16) provides us with a closed form expression for the A'"-electron average 
energy '^^'^^Eg,-y[nl^] as a function of the two-electron average energy '^^°'^^Eg,-y[nl'^] in 
the case of a general spin-independent two-body interaction. 

At this stage, it is convenient to use the (S, A)-parametrization defined by (2). In 
this parametrization, the energy '^^'^'^EQ_-,,[nl^] is a linear combination of the parameters 
E and A. The coefficient of S in this linear combination is clearly {1/2)N{N — 1). Since 
the parameter A appears only in the singlets of n£^, the coefficient of A in '^^'^^Eg,y[n£^] 
is (1/2) N{N — l)(nj/n/), where nj refers to the Young pattern [/] deduced from [/] 
by omiting its first row. (The latter assertion follows from the fact that Sq — for the 
singlet states so that S' = S and thus the sum on [/'] in (16) reduces to the orbital 
Young pattern [/] associated to the spin Young pattern [0\ = [ipi — 1, (fi2 — 1] for the 
spin S and the number of electrons N — 2.) Therefore, Eq. (16) may be written in the 
form 

- In{N- 1) + ^A^ . (17) 

The final step is to calculate the ratio nj/uf of the dimensions of the IRC's [/] (for 
Sn-2) and [/] (for Sn)- This may be achieved by using the well-known formula (see for 
instance Ref. [9]) for the IRC's of the symmetric group. We thus get 

n/ ^ ^2{^-i + 1) 

Uf N{N - 1) ■ ^ ^ 
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The introduction of (4) and (18) into (17) leads to Eq. (1). This completes the proof 
of the formula (1) for '^^'^^Eg.vi'ni^] in the general case of a spin-independent two-body 
Hamiltonian. If we are interested in the special case where V is the repulsive Coulomb 
interaction (ie, Vij = e^/rij), it is then sufl&cient to replace the parameters S and A 
in (1) by the appropriate linear combinations of the Slater-Condon-Shortley parameters 

3. Closing Remarks 

The main results of this paper [namely, Eqs. (16) and (1)] have been derived in 
the framework of atomic spectroscopy. Let us mention that the result (1) is also of 
relevance in the spectroscopy of partly-filled shell ions in condensed matter as far as S 
may be assumed to be a good quantum number (an assumption that is reasonable for 
transition-metal ions in crystals). The analysis leading to the results (1) and (16) rests 
on the use of the (single-configuration) shell-model and on the assumption of a unique 
single-particle radial wave-function for all the terms that the atomic (or crystal-field) 
configuration gives rise to. 

In the physical situation where 5 is a good quantum number (this situation occurs 
in some cases in atomic spectroscopy and in the spectroscopy of transition-metal ions 
in crystals), Eq. (1) provides us with an expression of the barycenter of the levels for 
a particular spin multiplicity ; then, the average energy for a given S corresponds to 
an observable and contact with experiment may be established. Another interest of the 
sum rule (1), valid even when S is not a good quantum number, is to be found in the 
fact that (1) is useful for the purpose of checking matrix elements when diagonalising 
the matrix of V by means of an electronic computer. As a third interest of (1), it is 
to be mentioned that the average energy (1) may be used in the method of spectral 
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distributions (mainly developed in nuclear physics) for the reconstruction of the global 
distribution of the energy levels with a fixed value of S. 

Let us observe that formulas which parallel Eqs. (16) and (1) might be derived 
in the shell model of nuclear physics. In this connection, we may think of obtaining 
average energies for fixed spin S, or fixed isospin T, or fixed spin S and isospin T. 
However, a complication arises when replacing N electrons on a n£-atomic shell by N 
nucleons on a n£-nuclear shell because the isospin degree of freedom manifests itself by 
the replacement of the trivial spin chain {U2)g D (81/2) g by the nontnma/ spin-isospin 
chain {U^)grp D {U2 — >■ SU2)g (8) {U2 — >■ SU2)rp (cf. Ref. [10]). In this respect, the use 
of class-sum operators [11, 12] might be of central interest. We hope to return on this 
extension to nuclear physics in a forthcoming paper. 

Finally, it is to be emphasized that the unusual sum rule (22) of the appendix is 
of pivotal importance in the derivation of Eqs. (16) and (1). This sum rule is based 
on Eq. (19) and on the Racah lemma. It is interesting to mention that Eq. (19) also 
plays a fundamental role in the derivation of sum rules for the intensity of two-photon 
transitions between Stark levels of a transition ion in a liquid or crystal environment 
[13]. 
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Let G be a finite or compact group. We use g to denote an IRC of G and to 
classify the rows and columns of a (standard) unitary matrix representation associ- 
ated to g. In this notation, the Clebsch-Gordan coefficients of the group G are written 
(^i^2A*iA*2bi^2^^A*) where b is an internal multiplicity label to be used when the Kro- 
necker product gi®g2 contains the IRC g several times. The Clebsch-Gordan coefficients 
of G satisfy ordinary unitarity relations controled by summations on /ii/i2 or bgjj. (see for 
example Refs. [14] and [15]). In addition, they satisfy the following sum rule (referred 
to as an orthogonality-completeness relation [15]) 



where /S.{g\gi® g2) is 1 or according to as g is contained in gi®g2 or not. In the special 
case where G is the group SU2 (or SO^)^ there is no need for the label h and we have g = j 
(or £) and /i = rrij (or mi) ; in this case, Eq. (19) is a simple rewriting of one of the two 
ordinary unitarity relations for the Clebsch-Gordan coefficients (jij2WiW2|iii2i"7,) that 
follows by using the symmetry property of the latter coefficients under the interchange 
h ^ j- 

Let us now consider a subgroup H of G. The label n may then be replaced by the 
triplet a/i7, where h stands for an IRC of the group H, 7 for an index to characterise the 
rows and columns of a (standard) unitary matrix representation Z)'* associated to h, and 
a for an external multiplicity label to be used when the IRC h of H occurs several times 
in the reduction of the IRC g of G. In the G D H basis, Eq. (19) may be transcribed as 

X] {gig2aihi'yia2h2'y2\gig2bgah'yy {gig2aihi'yia2h'2i'2\gig'2h' gah^i) = 




(19) 



A(^|^i®^2) 5{g'2,g2) S{ii2,H2) S{b',b) 



dim g 



dim g2 




H9\9i ® 92) S{g2, g2) 5(4, 02) 5(/t2, h2) (5(7^, 72) 5(6', h) 



dim g 




10 



The orthogonality-completeness relation (20) can be rewritten in terms of isoscalar fac- 
tors for the chain G D H. For this purpose, we use the Racah factorisation lemma 
[6] 

{gig20'ihi^ia2h2^2\gig2bgah^) = ^ {giaihi + g2a2h2\bgaph) (/ii/i27i72|^i^2/5/i7) 

(21) 

that gives an expression of a Clebsch-Gordan coefficient oi G, in a G D H basis, as 
a linear combination of Clebsch-Gordan coefficients (/ii/i27i72|/ii^2/3/i7) of H. The 
coefficients of this linear combination are the isoscalar factors (giaihi + g2a2h2\bga(3h) 
for the chain G D H. (In Eq. (21), the label /? is a multiplicity label of type b.) Then, 
by applying twice (21) in (20) and by using (19) for the group H in the so-obtained 
equation, we end up with 

^^^^ {giaihi + g2a2h2\bga(5h)* {giaihi + g2a'2h2\b' ga(5h) = 
^—f ^—^ '—f aimn.2 

a\h\ ah p (22) 

^{g\gi®g2) %2,^2) (^(4,02) 5{b\b) 

dim g2 

Equation (22) constitutes an orthogonality-completeness relation for the isoscalar factors 
of the chain G D H. 
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